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Abstract 

QCD predictions on the low- a; behaviour of the structure functions are reviewed and 
compared with the recent measurements of F2(x,Q 2 ) at HERA. The present theoretical 
accuracy of these predictions is discussed. 
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1 Introduction 



At present and future hadron colliders, precise quantitative tests of QCD and searches for new 
physics within the standard model (top and Higgs production) and beyond are carried out in 
the small-x kinematic region. By small x we mean that the ratio x = Q 2 /S, between the 
typical transferred momentum < p t >= Q in the process and the centre-of-mass energy yS 
of the colliding hadrons, is much smaller than unity. For these small- x processes reliable and 
accurate theoretical predictions are clearly necessary. 

Despite the general relevance of this topic, the present contribution does not intend to be 
a comprehensive (even if concise) review on the theory of small- x physics. I shall concentrate 
on the low- a; behaviour of the (proton) structure functions in deep inelastic electron-proton 
scattering. There are two very simple motivations for limiting ourselves to considering structure 
functions. First, the study is relevant by itself for investigating hadron physics at very high 
energy (since S = Q 2 /x, the small- x limit at fixed Q 2 is equivalent to S — > oo). Second, 
according to the parton model the structure functions are proportional to the parton densities 
of the hadron and their knowledge is essential for accurate phenomenology at any hadron 
collider. 

In the following, I shall review some of the main QCD predictions about the small-x be- 
haviour of the structure functions and, in particular, I shall try to discuss the present accuracy 
of these predictions. Some other aspects of small-x physics are considered elsewhere at this 
meeting [ED]. 



2 Structure functions at HERA 

The natural starting point of our discussion on low-x structure functions is represented by the 
recent HERA data on the proton structure function F 2 

HERA is the first electron-proton collider ever built and is operating at present with the 
centre-of-mass energy = 296 GeV. Due to this large centre-of-mass energy, it is able to see 
into very small values of the Bjorken variable x, x ^ 10~ 4 , even in the deep-inelastic-scattering 
(DIS) region 10 GeV 2 ^<5 2 ~10 5 GeV 2 (I use standard notations for the DIS kinematic vari- 
ables, namely Q 2 = —q 2 , x = Q 2 /2p ■ q, y = p ■ q/p ■ k with S = Q 2 /xy; see Fig. 1) This is an 
entirely new DIS regime which extends for two orders of magnitude, both in x and Q 2 , outside 
that (x^10~ 2 , Q 2 ^ 10 3 GeV 2 ) investigated by previous fixed-target experiments f5j. 

For values of Q 2 ^10 3 GeV 2 , the electron-proton inclusive cross section is dominated by the 
exchange of a single off-shell photon and given by 

F 2 (x,Q ) . (1) 



dxdQ 2 xQ 4 \ 2 l + R(x,Q 2 ) 



Here F 2 is the customary proton structure function, which is proportional to the sum <jl + 
(Jt, Ci's being the cross sections off a transversely (i = T) or longitudinally (i = L) polarized 
photon. The function R denotes the ratio R = olI&t'- it is smaller than 10% in the region which 
is being studied at HERA. It follows that the HERA data on the cross section ([!]) provides a 
direct measurement of F 2 (x, Q 2 ). 

The data collected at HERA in 1992 |pJQ] have already shown clearly that F 2 (x, Q 2 ) increases 
very steeply at small x. This behaviour has to be contrasted with that measured at higher re- 
values (Fig. 2). Previous data in the region of x ~ 10~ 2 are indeed consistent with a constant (or 
almost constant) behaviour of F 2 for x — > 0. Moreover, HERA data on F 2 as to be regarded as 
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a qualitatively new (although, possibly, not unexpected) experimental result: a similar strong 
rise with the energy has not been observed so far in total hadronic cross sections. 



Fig. 1: Kinematic variables for deep inelastic electron- proton scattering. 



Fig. 2: HI data on F 2 (x,Q 2 ). Data points of the NMC and BCDMS experiments are shown 
for comparison. 
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3 Total hadronic cross sections 

A compilation of data || on total hadronic cross sections is reported in Fig. 3. We can see that 
they increase slowly with increasing y/S. 



Fig. 3: Data for total cross sections with fits of type (0). 



The curves in Fig. 3 are the results of a Regge-type fit performed by Donnachie and Land- 
shoff fTj. It provides a very simple and economical description of all data in the form: 

a TOT = X S £ + Y S~ v . (2) 

The parameters X and Y depend on the type of colliding hadrons whilst the powers e and r\ 
are constrained to be 'universal' and found (by the fit) to be e = 0.08 and rj = 0.45. In this 
Regge-type phenomenological model «p(0) = 1 + e and 1 — 77 are thus respectively interpreted 
as the intercept of pomeron trajectory (it carries the vacuum quantum numbers) and of the 
p-meson trajectory (p, u>, f,a,). 
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The increase of the cross section (H) with the energy is due to the pomeron exchange and is 
indeed very slowly because the pomeron trajectory is above unity be a small amount, i.e. e = 
0.08. Actually, since e < 1, to a good approximation we have S £ ~ 1+e In S+ \e 2 In 2 S. This is 
the reason why double- logarithmic expressions of the type otot = AS~ n + B + C In S + D In 2 S 
|J are equally successful in describing the data. 

It is worth noting that this slow (approximately logarithmic) increase of the cross section 
holds true also for the photoproduction case. HERA data on the 7p cross section || (Fig. 4) 
are in remarkable agreement with the value predicted by Donnachie and Landshoff almost ten 
years ago |§. 



Fig. 4: Fit of type (§) for the 7p cross section data below v^S" = 20 GeV. The two data points 
at = 200 GeV correspond to the HI and ZEUS measurements at HERA. 



On the basis of this successful phenomenological description of the energy dependence of 
total cross sections, one can assume a corresponding x-dependence for the structure function 
F 2 HlO| , pA| . In particular, one can extrapolate the Regge-type parametrization @ from the 
photoproduction limit (Q 2 = 0) to the low-Q 2 {Q 2 ~ 10 GeV 2 ) regime by using the following 
expression [[II] 

F2{x ^ )=Ax ^^_ + B ^_^_. (3) 

This extrapolation is obtained by a simple dipole model (a ~ 750 MeV) for the Q 2 -dependence: 
F 2 is forced to vanish for Q 2 — > (gauge invariance) and to be Q 2 -independent at high Q 2 
(Bjorken scaling). 

The comparison between the parametrization (|3]) and HERA data (Fig. 5) shows that this 
type of models, based on the energy behaviour of low-p t hadronic interaction (soft processes), 
fails in describing the structure function F 2 (x, Q 2 ) as measured in the DIS regime. Hard physics 
is necessary for explaining the x-shape: a steeper F 2 is indeed expected from hard (perturbative) 
QCD. 
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Fig. 5: Comparison of the measured F2(x, Q 2 ) (HI data) with the parametrization in Eq. (|3|) 
(DOLA) and with several QCD expectations (see Sect. 5). The overall normalization uncer- 
tainty of 8% is not included in the error bars. 



4 Structure functions in perturbative QCD 
4.1 The QCD improved parton model 

According to the naive parton model the hadron is built up by point-like partons (quarks 
and gluons) which are almost non-interacting and scattered incoherently when their transverse 
momentum p±i is larger than some value Qo of the order of the inverse hadronic radius 1/Rh ~ 
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1 GeV. Therefore, in the large-Q 2 limit (Q 2 3> Ql) the structure function F 2 (x, Q 2 ) is given by 

f 2 (^ 2 )=Ee, 2 [/Jx)+/,(x)] , (4) 



i=i 



where f qi ( x ) = xqi(x) are the scale-invariant quark densities of the hadron and are the quark 
charges. 

Perturbative QCD modifies this simple picture. The large k±-g&p (Ql <C k± <C Qq) between 
the incoming proton and the scattered quark can be filled up by QCD radiation of massless 
partons. They have a logarithmic spectrum dk\jk\ in transverse momentum thus leading to 
corrections of the type 

^dk\ Q 2 



Ql k]_ Q 



that can be large (0:5 In Q 2 /Ql ~ 1 when Q 2 ^> Qq) even if the partons are weakly coupled 
(as(Q 2 ) <C 1). The large logarithmic contributions in Eq. (|) have to be computed to all orders 
in as and their resummation gives rise to parton densities evolving with Q 2 (scaling violations) 
in a predictable way. 

The QCD improved parton model is summarized by the following factorization formula ]T2| 

F 2 (x,Q 2 )= J2 f-PUx/z^Q'^QDUz^QD + Oil/Q 2 ). (6) 



Here f a (x,Ql) are the input (non-perturbative) parton densities (f 9 (x) = xg(x) is the gluon 
density) at the scale Qq, F 2a accounts for the QCD evolution from Ql to Q 2 and is computable 
in perturbation theory and the 0(1/ Q 2 ) term stands for the so-called higher-twist contributions, 
i.e. contributions vanishing as inverse powers of Q 2 in the high-Q 2 regime. 

Note that the Q 2 -dependence is completely factorized in Eq. (§) and thus predicted by 
the perturbative QCD component F 2a . On the contrary, the x-dependence of F 2 (x,Q 2 ) enters 
Eq. (§) through the convolution of F 2a and f a with respect to the longitudinal-momentum 
fraction z. As a result, the x-shape of F 2 is determined both by the input parton density and 
by its perturbative evolution. For instance, if we consider a model in which both the input and 
the evolution have a power dependence on x, i.e. f a (x) ~ x~ x and F 2a (x) ~ x~ a , from Eq. (^) 
we see that the steeper one wins as x — > : F 2 (x) ~ x -max{a,A}^ 

The perturbative QCD component F 2a in Eq. (§) is given by 

/•l J? 

F 2a (x, Q 2 ; Ql) = Y, - C»{x/z, *s(Q 2 )) F ba (z, Q 2 ; Q 2 ) , (7) 

b Jx Z 

where C 2a is the process-dependent coefficient function and F a b , the universal (i.e. process 
independent) Green function for the evolution of the parton densities, fulfils the following 
generalized Altarelli-Parisi equation 

r Q 2 Ah 2 pi 

F ab (x, Q 2 ; Ql) = 6 ab 5(1 - x) + £ / -± / cfe P ac (z, a^i)) F cb (x/z, Q 2 ) . (8) 

The parton densities at the hard scale Q 2 are obtained from the the corresponding inputs as 
follows 

/.(*, Q 2 ) = E T - G 2 ; «o) AC*. ^0) • (9) 

t z 
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The coefficient function C 2 (z,as) and the Altarelli-Parisi splitting function P(z,as) are 
both computable in QCD perturbation theory and their power series expansion has the form: 



C 2 (z,a s ) = 6(1 - z) + a s C?\z) + a 2 s C? } (z) + . . . , 
P(z, a s ) = a s P (1) (z) + a 2 P^ (z) + a 3 s P^ (z) + . . . , 



(10) 



(11) 



The coefficient function C 2 measures the effective coupling 7*-parton whereas Eq. (||) has 
a simple partonic interpretation. It describes a generalized parton evolution along a space- 
like cascade in which the successive fcj_'s are strongly ordered (Fig. 6). The Altarelli-Parisi 
splitting function P a b(z,as) can be consistently interpreted as the probability of emission of 
a bunch of partons (jets) whose transverse momenta are of the same order. In particular 
p(!) ; p( 2 ) ) . . . ; p( n ) ; in Eq. (|TI|) respectively represent the emission of one, two, ... , n partons 
with almost equal transverse momenta (Fig. 6b). We see that the radiation of an additional 
parton without /cj_-ordering costs an extra power of as and, hence, is a subdominant effect with 
respect to the Q 2 -evolution of F 2 (x,Q 2 ). 



Fig. 6: Partonic picture of the perturbative QCD evolution: (a) strong ordering of transverse 
momenta (Q 2 > kj_ n > . . . > k 2 ^ > Qq) along the space-like cascade; (6) splitting function for 
the emission of partons with comparable transverse momenta. 

For the sake of simplicity, in the following I do not consider the higher-order contributions 
( n ^ 1) to the coefficient function. Note, however, that by no means this implies that 
their quantitative effect is negligible in the small-x region [|T3| , |l4| , [l5| . 

In order to proceed in our discussion, it is also convenient to introduce iV-moments as follows 




(12) 



The limit x 



corresponds to the small- iV limit in the iV-moment space. 
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The Mellin transformation in Eq. fll2l) diagonalizes the convolution over z in Eq. @. The 
solution of the latter is 

f rQ 2 dk 2 1 
^(Q 2 ; Ql) = exp | M<*s{*i±)) \ . (13) 

where the anomalous dimensions 7 at are related to the iV-moments of the splitting functions 
as follows 

7a^(«s) = / ^ P(z, as) = a s Tn + a s iff + • • • • ( 14 ) 
Jo 

The expression ( |13D explicitly resums the large logarithmic contributions Og{hiQ 2 / 'Ql) k we 
have mentioned before. Inserting the one- loop (7$ ), two-loop (7^), •■• anomalous dimensions 
into Eq. (0), one can systematically resums leading (k = n), next-to-leading (k = n — 1), ... 
logarithmic terms. 



4.2 Altarelli-Parisi evolution to leading order 

The perturbative QCD evolution of the parton densities to leading order is particularly simple 
in the small-x limit |]T2j| . In this case the gluon channel dominates because the splitting function 

Pgg ( z ,a s )^^l}- + 0(a 2 s ), (s-0), (15) 

has a 1/z singularity, related to the exchange of a spin-one particle, the gluon, in the t-channel. 
The corresponding anomalous dimension (as = Cao.s/^) 

1uA<*s) * ^ + 0(a 2 ) , (JV - 0) , (16) 

is singular for iV — > 0. Inserting Eq. flTB"!) into Eq. (|13|) and transforming back to the x-space, 
one finds: 

F gg (x, Q 2 ; Ql) ~ exp |2 ^a 5 lniln^| . (17) 

The result in Eq. (|17|) shows that, due to the perturbative QCD evolution, the parton 
densities raise rapidly for x — > 0. Their increase is faster than any power of lnx, although 
slower than any power of 1/x. Note also that, because of the factor In Q 2 /Qq, the effective 
slope in x increases proportionally to the evolution range in k\. 

Equation (O) resums the double logarithmic contributions (as In - In %?) n to all orders in 

as. It has been obtained by assuming a fixed coupling as- The inclusion of the running 

(2) 

coupling and of the next-to-leading term 7^ for the anomalous dimensions does not change 
qualitatively the small-x behaviour of F(x, Q 2 ',Ql) (see Sect. 5). 

In order to gain physical insight into the result fll7|) , let us consider a simplified derivation. 
The probability of emission of a gluon with a large rapidity y = In 1/z and transverse momentum 
k± is: 

dk 2 

dw = as dy . (18) 

k ± 

The Green function F gg is obtained by integrating the spectrum ( Jl8| ) over the rapidity range 
L = In 1/x and summing over any number n of final state gluons as follows 

F gg (x,Q 2 ;Q 2 ) ~ E^ n sf L dy^---[ L dynf^---^Q(Q 2 >kl n >...>kl 1 >Ql) 

n Til JO JO J "'_Ln 

= H-^sLT [^...^e(Q 2 >k 2 ±n >...>kl 1 >Q 2 ) . (19) 
n n. j k ±1 K ±n 



s 



The statistical factor 1/n! in Eq. ([19|) is due to the identity of the n final-state gluons. The 0- 
function denotes the constraint of transverse-momentum ordering which is appropriate for the 
Q 2 -evolution to leading accuracy. After /^-integration, this constraint produces an additional 
factor of 1/n!, thus leading to ( (n!) 2 ~ (2n)!/2 2n for n ;§> 1) 

F 99 (x,Q 2 ;Ql)^Y,(^y (y/*sL HQVQl))^ - «p {2 Ja s L HQVQl)} • (20) 

This discussion shows that the rise of the parton densities at small x is due to the perturba- 
tive emission of many gluons (each of them giving a constant amplitude) over a large rapidity 
gap. Actually, at very small values of x, the L factors associated with the large rapidity gap 
may dominate the analogous factors \nQ 2 /Ql due to the fc^-evolution. In this case one has 
to go beyond the leading logarithmic approximation in Eqs. (|!5|)-(|l7l). The higher-order terms 
p( n \z) (n > 1) in the splitting functions ( [TT]) can be large and have to be taken into account. 



In fact, the additional powers of as due to the emission of gluons with comparable fcj_'s can 
be compensated by enhancing factors lnl/2;. Equivalently, we can say that we must give up 
/c_i_-ordering and compute the large In x-factors associated with the large rapidity gap and any 
value of transverse momenta. 

The possible impact of this improved treatment in the small- 2 region can be argued from 
Eqs. (|19D and (^0|) . Releasing the fcj_-ordering constraint in Eq. fll9|), we miss a contribution 
to Eq. (p0|), thus obtaining 

F gg (x, Q 2 ; Ql) ~ -r (const. a s L) n ~ exp {const. a s L} , (21) 

n n - 

where the constant factor in the exponent is related to the transverse-momentum integrations. 
As I shall discuss in the next subsection, the behaviour in Eq. (|2lD is precisely the result of the 
resummation of the terms (a^lnx)" - in perturbation theory. 

4.3 The BFKL equation 

In order to perform this resummation it is convenient to introduce the unintegrated parton 
density T{x, k; Qq) fTEfl . It gives the Green function for the QCD evolution with a fixed total 
transverse momentum k in the final state. The customary Green function F(x,Q 2 ;Qq) in 
Eqs. (0) and @ is then recovered by /c^-integration: 

F(x,Q 2 ;Ql)= [ Q d 2 k T(x,k;Q 2 ) . (22) 
Jo 

In terms of this function the resummation of leading contributions (as lnx) n is achieved by 
an integral equation whose kernel describes parton radiation without /cj_-ordering constraint: 

T(x, k; Ql) = TT 8(1 - x) <5(k 2 - Ql) + a s I' — f d 2 k' K(k; k') F(x/z, k'; Q 2 ) . (23) 

J x Z J 

The explicit form of the kernel is the following 

/ d 2 k! K(k- k!) T(x, k') = / ^4 k + q) - 0(k 2 - q 2 )^(x, k)l (24) 

J J nq* 1 J 



and was derived by Balitskii, Fadin, Kuraev and Lipatov (BFKL) almost two decades ago [17 . 
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The detailed discussion of Eq. (|23| ) within the context of the present paper can be found 
elsewhere [[Tj|[lI],[2(J . Here I limit myself to recall its main features. 

Since the BFKL equation (p3]) does not enforce /c^-ordering, it is not an evolution equation 



in kj_. Its general solution is thus contamined by higher- twist contributions. Nonetheless, in 
the large fc_|_-regime, one can disentangle the leading-twist behaviour in the form (Tn are the 
A^-moments of T(x)): 



7tv(«s) / k 2 



7jv(«s) 



Trk 2 



1 + 0(Ql/k 2 



(25) 



This power behaviour in kj_ implies that the resummation of the (as lnx)-terms performed 
by the BFKL equation is consistent with an effective gluon anomalous dimension r y g g t N( ( ^s) — 
7n(®s) |19ip0jpll ■ Its explicit expression is obtained by solving the following implicit equation 



as_ 
N 



x{lN{as)) 



(26) 



where the function xil) is related to the eigenvalues of the BKFL kernel and given by (ip(z) is 
the Euler -^-function) 

x(T)=2V(l)-V'(7)-V'(l-7) • (27) 
The power series solution of Eq. ( f2~6"| ) is (C(^) is the Riemann zeta function): 

o^^mf)"^ + 2C(3,(f)V 2 c (5 )(f) 6 + ((f)') . m 

The one-loop contribution on the r.h.s. of Eq. ( p8| ) reproduces the small- A" limit ( pTj[ ) of the 
Altarelli-Parisi splitting function P^ in Eq. (|ll|). The higher orders are the result of the 
resummation of the leading logarithmic contributions (agin™ -1 z)/z in P gg (z,as). 



Fig. 7: The BFKL characteristic function xil) f° r < 7 < 1. 7^ is the BFKL anomalous 
dimension. 



The characteristic function xil) m Eq. (p7|) has approximately a parabolic shape (Fig. 7). 
As x — > 0, A" decreases and reaches a minimum value N min = 4«sln2 at which 7^ has a 
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branch point singularity. Therefore the resummation of the singular terms (as/N) n builds up a 
stronger singularity at iV = iV m i n . This singularity, known as the perturbative QCD (or BFKL) 
pomeron, is responsible for a very steep behaviour of the QCD evolution at small x (and large 
k\): 



7Tk 2 \ Q 



x -4a s ln2 ? ^ 




F gg (x,Q 2 ;Q 2 ) = / d^H^Qi) ~ \ ^2 ^ S • (3°) 

I have so far reviewed the main qualitative predictions of perturbative QCD about the 
small- x behaviour of the parton densities and, hence, of the structure function F 2 (x, Q 2 ). The 
latter is expected to increase faster that any power of lnl/x. To leading logarithmic order in 
In Q 2 /Q 2 ) , the parton densities are driven by the Green function in Eq. (|TTD. Moreover, at very 
small values of x, the complete (i.e. including subleading terms in \b.Q 2 /Qq) resummation of 
the contributions (a^lnx)" - gives rise to the power increase in Eq. fl3"U|). Note that the transition 
from the small-x behaviour in Eq. (K7h to that in Eq. (R0) is also accompanied by strong scaling 



violations (compare exp JlnQ 2 /Q 2 ] with \/Q 2 /Q 2 ) 



4.4 Comment on unitarization 

Before presenting a more detailed quantitative discussion, I should add some comments on 
unitarity. 

Because of unitarity, the DIS cross section (|1|) cannot increase indefinitely. At asymptotic 
energies it must approach a constant (modulo In S'-corrections) limit as given by the geometrical 
size of the hadron. This limit implies the following unitarity bound on the structure function 
(Rh is the hadron radius) 

l^F 2 (x,Q 2 )<nR 2 h . (31) 

If we consider the leading-twist factorization formula @ and the fact that gluon evolution 
dominates at small x, from Eq. ([H]) we obtain a corresponding constraint on the gluon density, 
namely 

a s f g (x,Q 2 )<nR 2 h Q 2 . (32) 

On the other hand, from Eqs. (|17D and fl3~0|), we see that both the leading-order Altarelli-Parisi 
equation and the BFKL equation violate the constraint ( |3"2"D at sufficiently small values of x. 
This contradiction has a simple physical explanation. As soon as the gluon density becomes very 
large, the partons tend to overlap and are no longer scattered incoherently by the hard photon. 
In this regime we should take into account parton rescattering and parton recombination: they 
are responsible for fulfilling the unitarity bound in Eq. fl3~l|) . 

From a formal viewpoint, this means we have to include higher-twist contributions. Parton 
recombination is a process of order 1/(R 2 Q 2 ) and the probability of recombination of two 
gluons is approximately given by as{f g {x,Q 2 )) 2 . Therefore this process produces an higher- 
twist correction of relative order asfg^x^Q^/^R^Q 2 ) on the r.h.s. of Eq. (||). As soon as 
this corrections approaches unity, higher-twist terms cannot be neglected in the factorization 
formula @. We have thus recovered the constraint (p2[), which has to be correctly interpreted 
as a limit on the validity of the leading-twist approximation. 
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This argument shows that the perturbative QCD approach to small-x physics is (at least) 
self-consistent as regards to unitarity. The factorization formula (0) has to be used in a proper 
way by including higher-twist corrections whenever they are (estimated to be) large. 

In spite of many efforts and partial results jm , a general theory of rescattering has not 
yet been set up. Reasonable physical models are however available [fLSfl . Recent numerical 



estimates |2~3"] based on these models lead to the conclusion that rescattering effects give a 
small contribution (at most 10%) to F 2 in the HERA kinematic range. On the other side, at 
present, no evident signal of unitarity corrections (i.e. flattening of F 2 at small x) emerges from 
the HERA data. This is not surprising since, even for values of Q 2 as small as 1 GeV 2 , the 
bound fl31~l) gives F 2 ^25 (vri?| ~ 1 fm 2 ). Because of these reasons, in the following I shall not 
discuss any longer the issue of unitarity. 



5 NLO parton densities 

In order to perform a quantitative comparison between HERA data and perturbative QCD 
predictions, I shall consider some of next-to-leading order (NLO) parton densities that were 
available before HERA started operating. These parton densities are obtained by evolving (as 
in Eq. @) some input distributions according to the generalized Altarelli-Parisi equation (§) in 
NLO, i.e. including the full one- and two-loop [23| splitting functions (anomalous dimensions) in 



Eq. (|TT|) (Eq. (|14D ). The input parton densities at some scale Ql are fixed by fitting pre-HERA 
data on muon and neutrino DIS, prompt-photon production and Drell-Yan type processes 
in hadron-hadron collisions. The ensuing predictions have a high (better than 10%) overall 
accuracy, including the theoretical uncertainties due to higher-twist effects and perturbative 
corrections beyond two-loop order. For this reason, they provide a good starting point for 
detailed quantitative investigations of the small-x region. 

Actually, since pre-HERA data concern the high-x region, the input gluon and sea-quark 
distributions f g and / sea are essentially unconstrained for values of x smaller than 10~ 2 . The 
representative sets of NLO parton densities listed below precisely differ on the assumptions 
about the small-x behaviour of the input densities. 

The GRV parton distributions PS] are obtained by starting the perturbative QCD evolution 



at a very low input scale, Ql = 0.3 GeV 2 . At this scale the input distributions are chosen to 
vanish for x —>■ 0: f g {x) ~ x 2 , f se a.{x) ~ x ' 7 . Because of this 'valence-like' behaviour of the 
input, any steeper behaviour at high Q 2 is entirely generated by perturbative QCD evolution. 

The input scale used in the MRS analysis [^(| is Ql = 4 GeV 2 . They provide two rep- 
resentative sets of parton densities. Set D'o is obtained using input densities which are flat 
at small x: f g (x) ~ / SGa (^) ~ const. . Set D'_ instead uses very steep input distributions: 
f g (x) ~ fsea,(%) ~ x -0 ' 5 . Ro uglily speaking, the flat input corresponds to the so-called 'critical- 
pomeron' (with intercept equal to unity) behaviour whereas the steep input is reminiscent of 
the BFKL pomeron. 



As in the case of the MRS analysis, the CTEQ Collaboration |27| starts the perturba- 



tive evolution at the input scale Ql = 4 GeV 2 . Among the different sets of parton distri- 
butions they consider, set CTEQ1MS is that which has the steepest input densities, namely 
f g (x) ~ x _0 ' 38 (ln 1/x) 09 , f sea (x) ~ x -0 ' 27 . This small-x behaviour is just the result of a certain 
parametrization of the input and is not related to any 'hybrid pomeron'. 

Figures 5 and 8 show the comparison between the 1992 data| of the HI and ZEUS Collabo- 



^The comparison with the 1993 preliminary data of the ZEUS Collaboration, presented at this meeting [Q, 
does not substantially modify the discussion that follows. 
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rations and the structure function F2 as obtained from different sets of NLO parton densities. 
We see that the expectations from the GRV and MRSD'_ sets seem favoured by the data. 



Fig. 8: Comparison of the measured F 2 (x,Q 2 ) (ZEUS data) with the expectations from NLO 
parton densities. The overall normalization uncertainty of 7% is not included in the error bars. 

This observation, rather than being the main conclusion of the present contribution, gives 
me the opportunity for addressing some points which require further investigations and that I 
consider relevant for the understanding of small- x physics. 

The GRV parton distributions on one side and the MRSD' , MRSD'_, CTEQ1MS distri- 
butions on the other side represent, in a sense, two very different QCD expectations. 

In the GRV framework, due to the valence-like behaviour of the input densities, the fast 
increase of F2 at small-x and high Q 2 is entirely generated by perturbative QCD evolution as in 
Eq. (|17D . More precisely, according to Eq. (|17D , a structure function as steep as that reported in 
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Figs. 5 and 8 can be obtained only by means of a very long (Revolution. In the GRV case this 
is achieved by using a very low input scale. The weak (or, perhaps, amazing) point is that no 
higher-twist contribution is included (present) even at an input scale as low as Ql = 0.3 GeV 2 . 
Moreover, no small-x resummation (either from the BFKL equation or from higher orders in 
the Altarelli-Parisi equation) is considered. 

In the case of the MRS and CTEQ parton distributions the starting scale for the NLO QCD 
evolution is reasonably large, namely Qq = 4 GeV 2 . The different expectations for F<i at small 
x are entirely due to the different behaviour of the input densities. The MRSD'_ set, which is 
favoured by the data, has the steepest input densities. According to a widespread wisdom, the 
MRSD'_ framework is motivated by the BFKL equation. If one considers the BFKL pomeron 
intercept N min = 4asln2 ~ 2.65 as and a 'reasonable' value for as, say as — 0.2, one gets 
N m i n ~ 0.5: x~ ' 5 is indeed the small-x behaviour of the input parton densities in the MRSD'_ 
set. However, the perturbative (BFKL) pomeron is used here as non-perturbative input. On the 
contrary, the perturbative QCD evolution is performed with the NLO Altarelli-Parisi equation. 
In particular, since the input densities have an x-shape which is steeper than the one generated 
by the NLO QCD evolution, it follows (see the comment before Eq. (|7|)) that the small-x 
behaviour of F% is dominated by that of the non-perturbative input. I find this situation a little 
uncomfortable, although consistent at present. As a matter of fact, there are theoretical and 
phenomenological problems which are still to be solved. 

As regards to the theory, the BFKL equation resums the leading terms (as^x) n in QCD 
perturbation theory but it is not an evolution equation in k\. Its matching with the Altarelli- 
Parisi equation at larger values of x (or Q 2 ) and its extension (if any) to subleading orders 
in lnx requires the inclusion of the running coupling as(k 2 L ). Due to the lack of fc^-ordering, 
this leads to infrared instabilities coming from the Landau pole of as(kj_) in the low-fc^ (non- 
perturbative) region. 

As regards to phenomenology, in any perturbative expansion the leading term predicts just 
the order of magnitude of a given quantity. The accuracy of the prediction is instead controlled 
by the size of the higher-order corrections. 

The inclusion of subleading terms at small x is therefore necessary for precise quantitative 
studies. 



6 Beyond leading order at small x 



A possible approach [15,28| for improving and controlling the accuracy of the perturbative 
QCD predictions at small x (and large Q 2 ) consists in combining small-x resummation with 
leading-twist factorization. 

The generalized Altarelli-Parisi equation (§) for the evolution of the parton densities to 
leading-twist order systematically takes into account QCD corrections via the perturbative 
expansion (|TTJ) of the splitting functions P a b(z,as)- However, as discussed in Sects. 4.2 and 
4.3, higher-loop contributions to P a b are logarithmically enhanced at small x. More precisely, 



in the small-x limit, the n-loop order splitting function P^ behaves as 



P. 



1 r 



06 



x) ~ 



X 



\n n - L x + C(ln 



n-2 



Correspondingly, the anomalous dimensions 77V in Eq. (14]) have singularities for N 
form: 

'ocsy p (m 
N, 



lab,N(a S ) = £ (77) A aH + a S 
n=l L \ i V / 



n ab 



(33) 
in the 

(34) 
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These singularities may spoil the convergence of the perturbative expansion at small x. 
Nonetheless one can consider an improved perturbative expansion obtained by resumming the 
leading (A^), next-to-leading (B^), etc., coefficients in the small-x regime. Once these coef- 
ficients are known, they can unambiguously be supplemented with non-logarithmic (finite-x) 
contributions exactly calculable to any fixed-order in perturbation theory. Moreover, the effects 
of the running coupling can be consistently included according to Eq. ([13]), thus avoiding the 
infrared instabilities encountered in phenomenological attempts |f23| , |29|| to extend the BFKL 
equation beyond leading order. Obviously, in this approach, the absolute behaviour of the 
structure function F 2 at small x is still affected by that of the input parton densities. However, 
since logarithmic scaling violations are systematically under control, our confidence in making 
quantitative predictions at high Q 2 is enhanced. 

The present status of small- x resummation is the following. The gluon anomalous dimen- 
sions jgg^N contain leading singularities of the type (as/N) n . These contributions are given 
by the BFKL anomalous dimensions in Eq. (p8|). The next-to- leading corrections ots{cts/N) n 
in the gluon sector are still unknown. Beyond leading accuracy, the quark sector has to be 
considered on an equal footing with the gluon sector. The next-to-leading corrections to the 
quark anomalous dimensions have been computed recently ||28|| . The first perturbative terms 
are explicitly given by (Tr = 1/2 in QCD) 



%Mas) = ^ 5 .^ + f = ^p]|, (35) 




and higher- loop contributions can be found in Ref. [ 15 ]. 

A numerical analysis of the effect of next-to-leading corrections at small x has been presented 
recently by Ellis, Kunszt and Levin |3(J. They consider both the flat (D' ) and steep (D'_) input 



distributions used by MRS and, besides carrying out a full two-loop evolution (as in the MRS 



analysis), they include also the O ((a s/N) 4 ) -term of the gluon anomalous dimensions (28) and 
the O {{a s/N) 3 ) -term of the quark anomalous dimensions (0). Their results on the structure 
function F 2 are reported in Fig. 9. We can see that, for x ^ 10 -2 , the effect of the higher-order 
corrections strongly depends on the x-shape of the input distributions. The (Revolution of 
the steep input is marginally affected by the inclusion of higher-loop contributions^. On the 
contrary, the Q 2 -evolution of the flat input is pertubatively unstable^]. 

A first conclusion is that, if the input distribution is very steep, higher-order perturbative 
corrections are probably negligible, at least in the small-x region which is being investigated at 
HERA. On the other side, this does not mean that a flat input has to be excluded due to its 
perturbative instability. In this case one should rather include higher-loop corrections at small 
x and, possibly, recover reliable predictions after their all-order resummation. The difference 
between these two alternatives (steep input + fixed-order and flat input + resummation) is not 
just a matter of nominalism. Although the respective predictions for F 2 at a certain value of Q 2 
can be very similar, the scaling violations (i.e. the Q 2 - dependence) are stronger in the second 
case. Further studies are necessary for quantifying precisely the phenomenological consequences 
of this discussion for the HERA kinematic region, as well as, at future hadron colliders. 

*A similar perturbative stability is likely to occur in the GRV framework because of the very low input scale 
Qo and, hence, of the very long evolution in Q 2 . 

''Since in DIS the virtual photon couples directly to quarks (and not gluons), the quark anomalous dimensions 
(|35| ) (despite being formally subleading with respect to the gluon anomalous dimensions in Eq. (^8|)) are largely 
responsible for the perturbative instability in the HERA region. 
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Fig. 9: Results of the QCD evolution of the steep (upper curves) and flat (lower curves) 
input distributions. The effect of including higher-order terms (NLL) is compared with the full 
one-loop (L) and two-loop (NL) evolution. 

7 Conclusions 

The structure function F 2 (x, Q 2 ) measured at HERA increases very steeply for x<lCT 2 . A 
similar rise with the energy is not observed for total cross sections in soft hadronic processes. 
Therefore the structure function data at low-a; cannot be explained without hard QCD contri- 
butions. 

Perturbative QCD qualitatively accounts for the steep increase of F 2 . As regards to the 
sets of NLO parton densities available before HERA, the GRV and MRSD'_ parton densities 
are favoured by the HERA data. Both sets of parton densities are obtained by perturbative 
QCD evolution in two-loop order. In the GRV framework 'valence-like' input densities at a very 
low input scale undergo a very long Q 2 -evolution, thus acquiring the steep behaviour given by 
the Altarelli-Parisi equation QT7D . In the MRSD'_ framework, input distributions as steep as 
3.-0.5 (i nS pi re d by the BFKL equation) are considered at a much higher input scale. The low-x 
behaviour of F 2 at high Q 2 is thus dominated by that of the input densities^]. 

Estimating the present accuracy of the perturbative QCD predictions requires the evaluation 
of higher-order contributions and more detailed quantitative studies. For instance, from the 
discussion in Sect. 6, one can argue that flat (or almost flat, as suggested by the soft-pomeron 
behaviour 0) input densities at a scale of the order of 1 GeV 2 , combined with a dynamically 
generated perturbative pomeron (namely, all-order resummation in the Altarelli-Parisi splitting 
functions), may lead to structure functions consistent with HERA data. 

On the theoretical side, more calculations on next-to- leading terms at small x are necessary. 
On the phenomenological side, I consider relevant the analysis of the Q 2 -dependence and of the 
scaling violations. Studies on small-x effects in the structure of the hadronic final states are 

^The small- a; behaviour of the input densities can be 'tuned' to describing better F%. Recent fits to the 1992 
HERA data, carried out by MRS and the CTEQ Collaboration, give the behaviour a; -0 3 Q. 
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also important [0,0,^,0. More (and more accurate) experimental information will certainly 
be available in the very near future. 
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